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MODEL  APPLICATION 


The  model  selects  the  lowest  cost  combination  of  test 
equipment  configurations  and  test  periods  which  will  result 
in  the  desired  probability  that  the  item  tested  will  operate 
properly  when  used.  -It  is  particularly  applicable  to  items, 
such  as  missiles,  which  can  not  be  tried  to  see  if  they  work. 

1 The  model  will  be  most  useful  when  used  early  in  a program 

so  that  it  can  influence  test  equipment  design.  Nevertheless, 
if  the  test  equipment  has  already  been  purchased,  the  model 

can  be  used  to  determine  the  least  costly  test  periods, 

« 

The  model  can  be  used  to  simulate  up  to  three  echelons 
of  test.  These  might  be  a built  in  test,  an  intermediate 
test,  and  a depot  test.  The  model  automatically  indicates 
if  the  highest  echelon  of  test  is  unnecessary. 

The  model  can  be  used  to  evaluate  the  advantages  of  one 
location  of  a particular  test  echelon  over  another.  For 
example,  the  model  can  be  used  to  decide  which  is  the  best 
economic  choice,  shipboard  test  or  test  at  a weapons  station. 


.':cfaiN  m 1 

yr4 

tHk 
M t 

vtmmm 

■ 

a 

' JflnVNITNiU....— 

"TKr  >br.rFL>ga 

t IT  OlD 

j iimmiM/MiiuMim  «■  j 

Mt 

UULmt/mWrnUL  1 

\ 

•y . 

USING  THE  PROGRAM 


Input  Information 

To  determine  the  lowest  cost  test  scheme  one  must  start 
with  certain  information.  Table  1 lists  this  information 
It  has  been  divided  into  three  blocks.  Each  block  has  its 
own  input  format  which  will  be  discussed  below.  Some  explanation 
of  the  terms  may  be  of  value. 

The  cost  of  packing  and  shipping  to  the  ith  echelon 
tester  (Items  3 and  4 of  Block  1)  is  the  round  trip  cost  of 
packing  and  shipping  from  wherever  the  item  is  normally  used 
to  the  first  echelon  for  i=l,  from  first  to  second  echelon 
for  i=2,  and  from  second  to  third  for  i=3.  The  units  are 
dollars . 

The  phrase  "functions  tested"  (Items  7 and  8 of  Block  1) 
refers  to  the  tests  which  might  be  performed  on  the  item  to 
ascertain  if  it  is  in  good  operating  condition.  For  example, 
we  might  test  the  timing,  compression,  fuel  pres’sure,  battery 
charge,  spark  plug  gap,  and  breaker  point  gaps  of  an  internal 
combustion  engine.  Each  of  these  would  be  a function.  The 
cost  of  the  equipment  for  testing  the  compression,  for  example, 
would  be  the  cost  of  a compression  gauge.  The  time  required  to 
repair  this  equipment  would  be  the  time  needed  to  repair 
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the  compression  gauge.  The  failure  rate  of  the  equipment  to 
test  the  compression  would  be  the  failure  rate  of  compression 
gauges.  The  time  required  to  test  the  compression  would  be  the 
time  to  unscrew  the  spark  plugs,  insert  the  gauge  and  read  the 
pressure,  and  reinstall  the  plugs.  The  failure  rate  of  the 
compression  would  be  the  number  of  compression  failures  per  use 
cycle.  We  might  choose  one  hour  as  a use  cycle  or  we  might 
pick  some  more  convenient  period  depending  on  the  application 
of  the  engine. 

These  functions  must  be  ordered  according  to  increasing 
values  of  the  ratio  of  the  cost  of  the  equipment  needed  to 
test  the  function  to  the  failure  rate  of  the  function  tested. 
This  will  result  in  detecting  the  greatest  number  of  failures 
per  test  set  dollar  at  first  echelon.  Usually  there  are  more 
first  echelon  testers  than  second  echelon  and  third  echelon 
testers.  Therefore,  this  procedure  will  result  in  least  test 
equipment  cost.  The  order  number  of  a function  is  referred 
to  as  m. 

In  many  cases  one  will  have  some  idea  what  tests  should 
be  performed  at  each  echelon.  If  so,  limits  can  be  set  on  the 
range  of  the  m's  which  define  what  testS  are  performed.  The 
variation  of  total  costs  as  the  m's  change  within  the  range 
can  be  studied.  The  initial  and  final  values  of  the  index, 
m,  of  functions  tested  and  also  the  increment  by  which  the 
index  is  to  be  increased  each  loop  are  specified  under  Items 
7 and  8 of  Block  1. 
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The  number  of  functions  tested  (Item  1 of  Block  2)  is 
the  total  number  of  functions  of  the  item  which  it  is  proposed 
to  test.  It  is  possible  that  the  analysis  will  show  that  it 
is  not  economically  desirable  to  test  all  of  these.  The  value 
of  item  1 of  Block  2 must  not  exceed  20. 

The  fraction  of  items  in  use  which  are  not  defective 
(Item  2 of  Block  2)  is  the  fraction  of  the  time  that  when  an 
item  is  used  it  will  operate  properly.  This  is  sometimes  called 
the  operational  ready  rate.  It  is  the  major  constraint  in 
minimizing  testing  costs. 

The  repair  material  rate  (Item  3 of  Block  2)  is  the  ratio 
of  the  average  cost  of  the  materials  to  make  one  repair,  to  the 
item  cost. 

The  other  input  quantities  are  self  explanatory.  To 
simulate  two  test  echelons,  the  first  two  values  of  Item  8, 

Block  1 are  set  equal  to  the  total  number  of  functions  which 
can  be  tested  and  the  third  value  is  set  equal  to  one.  The 
number  of  testers  at  third  echelon  is  set  equal  to  zero. 
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INPUT  DATA  BLOCKS 
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of  functions  tested  (m<20)  (dimensionless) 


0) 

s 

3 

Z 


X) 

3 

-3 

X 

C 

3 

-3 

X 

3 

3 

3 

£^ 

X 

*3 

E 

-3 

> 

3 

-3 

X 

3 

O 

O 3 

3 

•H 

O 

O P 

r— 4 

P o 

> 

•3 

in 

t/i 

jC  3 

O 

W 

r— 4 

z 

O 

3 

in 

(/) 

4-4  O 

3 

3 3 

O 

t-4 

1—4 

s 

a> 

O 

f— 1 

3 

•H  P 

44 

m 

1—4 

o 

tH 

r^ 

3 3 

£ 3 

O 

3 

C 

c 

3 X >s 

O 

3 O 

£ 

3 3 

44 

3 

o 

o 

3.  O X) 

CM 

O -H 

O 3 

O 

•H 

•H 

3 

•H  Cfl 

o 

44  iH 

£ 

3 

in 

(/) 

••  / — \ 

P 3 

p 

•H 

C 

c 

tOTd  CM 

•• 

CJ  X 

3 3 

o 

£ 

a> 

0) 

3 3 3 

T3 

3 

f— 4 

3 3 

p 

e 

B 

O 

O 

3 3 

-3  ' 

o 

•1-4 

•H 

•>  O V4 

X 

<P  -H 

E 

3 

1—4 

p 

•3 

'd 

cn  3 O 

p 

3 

o 

O • 

■♦J  (fl  w 

3 

X 3. 

3 

P 

£ ‘ — > 

3 

in 

£ 

P 3 

44 

C 3 

O 3 

U 

3 H to 

•H  3 

•H  3 

3 3 

£ 

> 

l-^ 

4-»  3 3 

3 

3 

£ 

44  O 

O 

•H 

O- 

o 

P O 

3 

C 3 

Z 

3 

'd 

3 oo 

3 

in  P 

-3 

O -H 

3 \ 

3 

o 

d 

O - 3 

X 

3 

3 

•H  3 

3 m 

0> 

d 

1—4  \t4 

a. 

P -3 

o 

P O' 

-3  3 

3 

M-l 

3 CM3 

3 

3 

3 3 

3 3 

3 

(U 

o 

•r;  3 -H 

(3 

O 3 

C 

3 

O 3 

-3 

TJ 

/ — » 

0^6 

1 

P -H 

•H 

3 P 

r-4 

3 

c 

»-i 

3 ~ 3 /-V 

3 

3 

/— \ 

4-1  W 

C -H 

O 

>1-1 

0) 

cd 

in  ID  in 

o 

-3  O' 

C<o> 

3 

•H  3 

o 

V 

"3  >4-)  .(J  t/) 

p 

3 3 

O 

X P 

3 

C 

c 

> 

>- 

3 <n  3 3 

S ■ 

in  3 

•H 

P 

£'^ 

•H 

•M 

O 3 -3  1-4 

3 

3 

P 

•H  O 

3 

0) 

O 

in 

O >3  3 

z 

3 

o • 

•H 

P • 

3 

i-< 

•o 

at 

3 3 0 

P £ 

C p 

p m 

•H  3 

O 

(4) 

f— * 

cn  3 O -H 

3 

3 -H 

3 C 

m 3 

3 

•H 

C 

tj 

O Mh  w 

O 

3 P 

44  3 

3 X 

3 O 

P 

JS 

o 

X 

m tH  3 

P 

£ 

• E 

P 

O X 

O 

o 

•H 

u 

0 3-33 

P,  (U 

X 3 

44 

? 

3 

•H 

•M 

v— / 

-3  3 6 

(rt 

•H  X 

P -H 

O O 

c 

3 

rC 

U 

3 u cn  .i-c 

C 

3 P 

•H  3 

P 

O 3 

3 

? 

d 

3 3 3-3 

o 

O' 

cr 

3 

•H 

Cd 

X)  ^ 

•H 

3 

P 3 

P P 

p m 

X 

3 

0) 

6 4-)  w 

X 

m 

C 3 

u 3 

P 

V) 

D cn  -H  -3 

d 

3 

3 P 

3 3 

3 

•H  • 

3 

3 3 O 

•H  -O 

P m 

£ 

3 £ 

3 

o 

»-> 

•r-l  f-i  Xt 

o 

3 3 

3 

3 3 

3 3 

P £ 

C 

•H 

4> 

3 ‘4H  3 -P 

•p 

3,  > 

O P 

•H  3 

P 

in  'H 

•H 

P< 

X XI  3 

•H 

3 O (A 

P 

3 3 

X -H 

3 P 

d 

p 3 g 6 

3 t-H  3 

o 

O'.H 

P 

p 

tn 

in 

O 3 

<4^ 

tH  3 

P -H 

3 'H 

•H  3 

Pi 

£ 

V 

■p  ^ C 3 

O (/) 

O O O 

3 w 

3 

O 

O 3 

0) 

rH 

o o 

(/) 

P 44  Z 

3 3 

44  44 

3 

P 

■M 

d 

o 

P Cfl  cfl 

^ a> 

N— / 

£ X 

O 

O 3 

P 

•H 

•H 

3 en  -H  X 

o ^ 

-3  £ 

3 

X 

P 

-3  3 

o 

3 3 X P< 

X 3 

3 

•H  44 

3 X 

3 3 

3 in 

<4-4 

0) 

1— 4 P E“*  rt 

£ O 

3 O P 

3 O 

P 

P 3 

3 

O 

0) 

3 cx^ 

3 -H 

•H  P C 

O' 

3 -3 

3 

•H  P 

d 

in 

> C ' 

3 « 

3 3 

3 P 

3 3 

3 3 

3 in 

3 

€ 

3 

O • 3 

3 

O-rH  g 

in 

3 

o'  3 

o 

f— 4 f— 4 O 

£ 3 

3 3 

44  O 

3 O 

3 3 

3 P 

•H 

M 

3 3 P P 

3 g 

3 g -H 

o o 

3 1—4 

3 3 

3 

♦•4 

•H 

(U 

•H  X W ? 

£ -3 

O 3 

3 r-4 

3 -H 

3 

O 

d 

•Q 

P U 3 3 

•H  -3 

3 3 O' 

P X 

r-4  O 

3 3 

3 X 

cd 

a 

6 

•H  3 P Z 

X ^ 

£44  3 

l/tx 

•H  44 

•H  3 

£ P 

0) 

3 

C 

3 

•H 

o 

3 

3 

•H 

o. 

a: 

Z 

1— 1 

2 

H 

u 

3 

3 

H 

• 

• 

• 

• 

• 

• 

• 

K> 

m 

iH 

r4 

to 

ifi 

CM 

to 

X 

u 

i< 

o 

u 

X 

o 

CQ 

X 

n 

6 


Input  Format 


The  input  data  format  is  different  for  each  of  the 
three  blocks.  The  first  is  an  11  x 3 matrix,  the  second 
a 6 X 1 matrix,  the  third  a 5 x the  total  number  of  missile 
functions  being  tested  plus  one. 

BLOCK  1 - consists  of  11  cards  of  data,  3 items  per  card. 


Data  are  punched  in  successive  10  column  fields 
beginning  in  column  21.  The  format  of  a block  1 
data  card  is  given  below. 


uflnoBflBpngmisi^iaiiiMaasags^afflgadias^ 

sasa 

^9l6^ll6^j64i65jfi6j6T!6«W70!ol.7p’?*^:«'7V«  ’S®' 

fill 

tmi 

iiniiiiiiiiiiiliiiiiii 

BLOCK  2 - consists  of  six  cards  consisting  of  one  data  item 
each,  punched  in  a 10  column  field  beginning  in 
column  21.  The  format  of  a block  2 data  card  is 
given  below. 
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BLOCK  3 - consists  of  5 rows  of  data  each  of  which  contains 
m items  (m521) . The  m data  items  are  punched  4 
items  per  card  in  successive  15  column  fields  begin- 
ning in  column  21.  The  format  of  a block  3 data  card 
is  given  below. 


Since  m may  equal  21,  up  to  «five  cards  of  data  may  be 
required  for  each  row  of  block  3 data»  When  punching 
data  cards  always  use  the  decimal  points  in  each  data 
item. 

Figure  1 shows  how  the  deck  is  arranged.  In  the  figure 
a delta  indicates  blank.  The  punching  for  all  cards  except 
the  data  deck  begins  in  column  1.  The  region  value  on  the 
second  card  may  be  changed  to  accommodate  a large  number  of 
cases  of  n^^,  n2 , m^ , and  m2. 
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Interpreting  Results 


Table  3 sho\'s  the  results  of  a sample  run.  The  far 
right  hand  column  in  Test  Cost  Summary  is  the  cost  of  testing 
using  the  equipment  indicated  by  the  numbers  in  the  first  two 
columns  and  the  test  schedule  from  the  next  three.  An  I in  the 
fifth  column  indicates  that  the  third  echelon  test  was  not 
necessary  to  obtain  the  required  probability,  U,  that  the  item 
is  unfailed  when  needed.  A star  in  the  fifth  column  indicates 
that  for  some  reason  the  program  was  unable  to  calculate  a 
value  for  n3.  The  reason  is  also  given.  The  most  common 
cause  of  an  inability  to  calculate  is  that  the  third  echelon 
test  must  be  performed  more  frequently  than  the  second  echelon 
test  in  order  to  obtain  the  required  value  of  U.  Frequently, 
this  condition  can  be  detected  by  looking  at  the  preceeding 
values  of  n^  for  the  same  values  of  m^ , m2,  and  n-j^ . If  n^  was 
getting  close  to  one,  it  is  likely  that  it  would  have  to  be 
less  than  one  and  there  would  not  be  a solution.  In  either 
of  these  two  cases  no  action  is  necessary. 

The  sixth,  seventh,  and  eighth  columns  are  the  probabilities 
that  failed  items  will  be  detected  at,  or  before,  the  test 
level.  The  ninth,  tenth,  and  eleventh  columns  are  the  cost  of 
performing  a single  test  at  the  test  level,  or  echelon.  The 
twelfth  column  is  the  cost  per  use  cycle  of  testing  all  of  the 
items.  It  does  not  include  the  cost  of  the  test  equipment. 

Occasionally,  a solution  is  not  found  on  the  first  try 
and  it  is  not  obvious  that  the  reason  is  either  of  those  given 
above.  Possible  reaction  to  the  error  messages  are: 
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1.  NO  CONVERGENCE  WITHIN  ITERATIONS.  A solution 
may  not  be  obtained  because  more  iterations  are  needed 
to  produce  the  desired  accuracy.  The  situation  can 

be  corrected  by  increasing  the  input  quantity  "maximum 
number  o£  iterations  for  Newton-Raphson  method". 

2.  AT  SOME  ITERATION  STEP  DERIVATIVE  IS  ZERO.  If  a 
zero  is  encountered  in  the  denominator  of  an  expression^ 
the  Newton-Raphson  method  terminates. 

3.  INITIAL  X VALUE  WILL  CAUSE  EXPONENT  OVERFLOW.  NEED 
NEW  STARTING  VALUE.  Because  the  expression  being 
solved  contains  exponentials  it  is  possible  that  a 
value  too  large  to  be  handled  by  the  computer  will 

be  encountered.  This  difficulty  can  be  eliminated 
be  selecting  another  initial  value  for  Item  5 of  Block  2. 

The  range  of  should  be  expanded  if  the  service 
life  cost  does  not  pass  through  a minimum  for  every  combination 
of  m , m. , and  n . An  exception  to  this  rule  is  if  the  third 

^ ^ X 

echelon  test  would  have  to  be  performed  more  frequently  than 
the  second  echelon  test  as  discussed  above. 

Once  there  is  a minimum  service  life  cost  of  testing 
for  every  combination  of  m^^,  m2 » and  n^^ ; check  to  see  if  there  is 
also  a minimum  of  these  minimums  for  every  combination  of 
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nij^  and  ni2 . If  there  is  not,  increase  the  range  of  n^^  until 
a minimum  is  obtained.  The  minimum  of  this  minimum  of  the 
minimums  is  the  least  cost  testing  scheme.  That  is,  the  first 
m^  functions  should  be  tested  at  the  first  echelon,  the  next 
m2-mj^  functions  should  be  tested  at  second  echelon,  the  remaining 
functions  should  be  tested  at  third  echelon,  first  echelon 
tests  should  be  performed  every  n^  times  the  item  is  used, 
second  echelon  tests  should  be  performed  every  n^  times  a 
first  echelon  test  is  performed,  and  third  echelon  tests  should 
be  performed  every  nj  times  the  second  echelon  test  is  performed. 
Naturally,  third  echelon  tests  will  not  be  necessary  if  an 
I appears  under  nj. 

Table  3 is  a sample  run.  It  will  be  seen  that  as  n2 
is  varied,  holding  m^^ , m2»  and  n^  constant,  a minimum  was 
not  always  reached.  For  example,  when  m^^  is  1 , m2  is  2, 
and  n^  is  1 the  total  cost  is  still  decreasing.  One  might 
guess  that  it  will  never  get  as  low  as  the  $14,356,321  which 
is  reached  for  mi=4,  m2=10,  and  ^^^=5  but,  to  be  sure,  the  range 
of  n2  should  be  increased  until  a minimum  is  found. 

In  the  case  of  mj^  = l,  ni2  = 2,  and  nj^  = 5,  n3  does  become  less 
than  1.  In  cases  like  this,  the  lowest  cost  scherie  for  which 
nj  is  not  less  than  1 will  be  the  minimum. 
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In  the  case  of  mj^=2,  m2=10,  and  '^3"!  indicating 

that  third  echelon  testing  is  not  necessary.  The  service  life 
cost  is  calculated  assuming  no  third  echelon  test  equip- 
ment is  purchased.  For  these  cases,  the  range  of  n^  must  be 
expanded  until  a minimum  is  reached.  The  case  mj^=2,  m2  = 10, 
and  n =3  illiistrates  what  may  happen.  In  this  case  n,  becomes 
finite,  and  then  becomes  less  than  1.  It  turned  out,  in  this 
case,  that  a minimum  was  reached  first. 

We  must  also  be  able  to  find  a minimum  of  the  minimums  for 
every  set  of  m^^  and  m2.  In  at  least  two  cases,  not  only  is 
this  minimum  of  minimums  not  found  within  the  range  of  n^,  but 
it  seems  very  likely  that  the  minimum  cost  may  be  less  than 
the  lowest  cost  found  in  this  run.  These  cases  are  m^=l  and  m2=4 
and  6.  There  may  be  others.  Consequently,  another  run,  with 
increased  n^^  range,  is  necessary.  This  second  run  is  given 
in  table  4.  It  will  be  seen  that  the  range  of  n^^  is  still  not 
great  enough.  Therefore,  the  run  of  table  5 was  made. 

In  preparing  table  5 , the  indices  were  incremented  by  1 
because  it  was  expected  that  this  would  be  the  final  run  and 
it  was  necessary  to  consider  all  possible  test  equipment 
configurations  and  test  intervals .The  range  of  n^  was  extended 
to  complement  the  previous  results.  The  range  of  n2  was  contracted 
based  on  the  results  of  table  4. 


In  every  case  but  mj^  = 10,  in2  = ll  a minimum  is  found  or  at 
least  it  is  obvious  that  the  values  of  the  service  life  costs 
have  leveled  off.  For  the  exception  it  can  be  argued  that  this 


case  could  not  be  the  least  cost  since  the  values  have  been 
increasing  with  m^^ . It  will  be  seen  that  the  cheapest  combin- 
ation of  test  equipment  configurations  and  test  periods  is  for 
testing  every  14  use  cycles  at  first  echelon  with  a tester 
which  tests  the  first  5 functions,  testing  every  14X3=42 
use  cycles  with  a tester  which  tests  the  next  11-5=6  functions, 
and  testing  every  14X3X6=252  use  cycles  with  a tester  that 
test  the  last  12-11  = 1 function.-  We  have  not  considered  the 
possibility  of  performing  all  tests  at  first  and  second  echelon 
and  none  at  third.  It  is  possible  that  a lower  cost  might  be 
found  if  this  were  done. 

Table  6 is  a run  for  two  echelons  of  test.  It  will  be 
seen  that  a minimum  as  a function  of  n^^  is  not  found  in  this  run 
and  it  was  necessary  to  increase  the  range  of  n^^  as  shown 
in  table  7.  In  this  table  it  is  seen  that  the  minimum  cost 
of  testing  using  only  two  echelons  is  greater  than  for  using 
three  echelons.  Therefore,  the  minimum  cost  test  equipment 
configurations  and  test  periods  are  those  cited  above. 
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MATHEMATICAL  MODEL 


The  mathematical  model  is  based  on  the  testing  scheme 
which  is  diagramed  in  figure  2*  The  first  echelon  test  is 
performed  every  n^^  times  the  item  is  used.  The  second  echelon 
test  is  performed  every  n2  times  the  first  echelon  test  is  performed. 

The  third  echelon  test  is  performed  every  n^  times  the  second 
echelon  test  is  performed.  When  an  item  fails  a test  it  is 
repaired  or  replaced  and  introduced  as  a new  item.  After  the 
testing  scheme  has  been  in  operation  for  some  time  an  equilibrium 
condition  will  obtain  in  which  the  probability  that  an  item 
which  is  selected  at  random  will  be  unfailed  is  constant. 

We  wish  to  know  the  relation  between  this  constant  and  the 
parameters  of  the  test  arrangements. 

Since  we  start  each  use  cycle  with  some  number,  M,  of  good 
units,  the  number  of  good  units  which  have  been  used  x times  is 
given  by  the  expression 

-a.x 

U.Me 


where  is  the  fraction  of  new  units  which  are  good,  M is 
the  number  of  new  units  introduced  each  use  cycle,  and  a 
is  the  failure  rate  per  use  cycle  of  the  item.  At  the  beginning 
of  each  use  cycle,  some  units  have  not  been  used,  some  have 
been  used  once,  some  twice,  and  so  on  up  to  (nj^n2nj-l).  Items 
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that  have  been  used  n2^n2n2  times  are  removed  from  service  and 
tested  on  the  third  echelon  tester.  If  they  pass  this  test 
they  will  be  returned  to  service.  If  they  do  not  pass  this 
test  they  will  be  replaced  with  units  which  are  new  or  reworked. 
The  new  or  reworked  units  will  have  been  tested  on  the  third 
echelon  tester.  If  the  probability  that  a unit  which  passes 
this  test  will  really  be  good  is  Up  then  the  number  of  units 
which  are  in  service  which  are  good  will  be  given  by 

-eoL 

= E.Q1 

L—  I 

0 

Now,  M is  the  number  of  new  units,  or  units  which  have  passed 
the  third  echelon  test,  which  are  introduced  each  use  cycle. 

It  is  the  sum  of  the  number  of  units  which  fail  the  first  and 
second  echelon  tests  each  use  cycle  plus  the  number  of  units 
which  are  tested  on  the  third  echelon  test  each  use  cycle. 

The  number  of  units  tested  on  the  first  echelon  tester  each 
use  cycle  is 

Jl 

where  N is  the  total  number  of  units. 
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Of  these 


C- 


fail  the  test.  is  the  fraction  of  failed  items  which  pass  the 

first  echelon  test.  The  number  of  units  which  survive  to  be  tested 
on  the  second  echelon  test  the  first  time  is 

na 


Cf  these 


fail  the  test.  FA2  is  the  fraction  of  failed  items  which  pass  the 
second  echelon  test.  The  number  which  do  not  fail  will  be 

lit 
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Of  these 


survive  to  be  tested  on  the  second  echelon  tester  the  second 
time.  Of  these 


r\<L. 


-CX-. 


fail  and 

n <2, 


survive.  Extending  this  reasoning,  one  finds  that  the  total 
number  of  items  rejected  each  use  cycle  by  the  second  echelon 
test  is 

T -v^FA, 7 

M <2.  Z-  «-  L'  j 
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The  number  of  items  which  survive  to  be  tested  on  the  third 
echelon  tester  is 


-a.(\-p^z)  n.n^hj 

H e. 


Therefore,  the  number  of  new  units  introduced  each  use  cycle 
will  be 


M r -a/i-P^i)n.|  . 


+ I1e 


-o^i-FAi)K],ni  hj 


l-(L 
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By  substitution  we  get  the  fraction  of  the  items  which  are 
good  at  any  time: 


U- 
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u.  1 

i-e.  J 

.1-  e. 
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- 

1- 

a 
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Eq.2. 


Normally,  a test  scheme  will  be  designed  to  provide  a U 
which  is  greater  than  or  equal  to  some  specified  value. 
Obtaining  greater  values  of  U increases  costs.  Decreasing 
the  false  acceptance  rates  increases  U and  adds  to  the  cost 
of  the  test  equipment.  Decreasing  the  values  of  the  n's 
increases  U and  increases  labor,  shipping,  and  spares  costs. 

There  is  a selection  of  false  acceptance  rates  and  of  test 
intervals,  n's,  which  provides  the  required  U at  minimum  cost. 
We  must  find  this  selection. 

The  number  of  tests  which  will  be  performed  at  each  echelon 

each  use  cycle  depends  on  the  values  of  n^^ , n , and  n and  on 

2 3 

the  false  acceptance  rates.  The  number  of  items  which  will  be 
tested  at  the  first  echelon  each  use  cycle  is 

h, 

The  number  of  units  which  will  be  tested  at  the  second  echelon 
for  the  first  time  is 
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. , -o.(i-FAi)  n.n, 

M e. 


and  for  the  second  time  is 

. -0uti-FA,)ah,nx  -o-(FA.-'=«Oi^-n^ 
Me.  e. 


The  total  number  of  units  tested  at  the  second  echelon  will  be 

-ft.(«-FA,)  riiin^  I— (e.  

I --e. 


each  use  cycle.  Now  since 
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this  becomes 


N 


-<X.(|-Frt, 


r_ 


1 


The  number  o£  items  tested  each  use  cycle  at  the  third  echelon 


will  be 


n e. 
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1 

[l-e  j 

1 

1-  e. 

Now,  let  be  the  costs  associated  with  performing  tests  on 
each  unit  at  the  first  echelon,  K2  the  cost  at  the  second  echelon, 
and  the  cost  at  the  third  echelon.  The  cost,  K,  of  testing 
incurred  each  use  cycle  will  be 


1 ' 
' 'ri. 


N ^ 
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This  quantity  times  the  total  number  of  use  cycles  in  the 
service  life,  when  added  to  the  cost  of  the  test  equipment,  will 
be  the  total  service  life  cost. 

The  cost  of  a test  set  to  achieve  the  required  false 
acceptance  rate  is  determined  as  follows.  A test  set  is  required 
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which  will  detect  all  but  FA  o£  the  failed  items  tested. 

We  can  choose  the  functions  to  be  tested  to  minimize  the 
cost.  Let  a^  be  the  failure  rate  of  the  ith  function  which 
the  item  performs.  Let  be  the  cost  of  the  equipment  to 
test  each  of  the  m functions. 


If  there  are  a total  of  mj  functions  which  might  be 
tested  and  if  it  costs  for  the  equipment  to  detect  the  failure 
of  the  mth  function,  the  cost  per  fraction  of  failures  detected 
will  be 


- rA 


C Vv, 


where  is  the  fraction  of  the  failed  items 
the  test  if  the  mth  function  is  not  tested  and 


which  will  pass 
is  given  by 


■*  - ' 

where  a is  the  failure  rate  of  the  item. 


of 


If  the  functions  are  ordered  according  to  increasing  value 


(X, 
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they  will  be  in  the  order  in  which  they  should  be  chosen  to 

« 

obtain  the  required  test  set.  The  false  acceptance  rates  of 
the  first  and  second  echelon  testers  are  given  by 

•VM 


FA,*  I - 


I- 


where  is  the  highest  index  of  the  functions  tested  at  first 
echelon  and  m2  is  the  highest  iijdex  for  second  echelon. 

The  false  acceptance  rate  of  the  third  echelon  tester  is 

l-Uo-  1 - 

Ou 

The  cost  of  the  first  echelon  tester  will  be 

rv\| 


the  cost  of  the  second  echelon  tester  will  be 
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and  the  cost  of  the  third  echelon  tester  will  be 

, te..  C,  , = C -t-  ) 

where  is  the  cost  of  the  basic  tester. 


The  cost  of  performing  each  test  is  the  cost  of  the 
labor,  the  shipping  to  the  test  site,  packing  for  shipping, 
the  spares  to  replace  the  items  while  they  are  away  for  test, 
the  test  set  maintenance,  and  the  consumables  associated  with 
the  test.  These  costs  are  given  by 


'h^ <x.r-£>v  ^ d 1 Vice.  Li 


fid 


where  the  subscript,  i,  indicates  the  echelon  of  the  test. 


Some  of  the  quantities  which  depend  on  the  echelon  of 
test  depend  on  what  tests  are  performed  and,  through  this 
dependence,  on  the  false  acceptance  rate.  These  quantities 
are  the  test  time,  the  test  sets  MTBR  and  MTBF,  and  the  average 
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materials  cost  for  repairing  the  test  set.  These  quantities 

can  be  related  to  the  functions  which  the  test  set  tests. 

For  example,  it  will  take  a particular  time  to  test  a particular 

function.  Consequently,  the  quantities  can  be  found  as  follows. 

The  test  time  is  the  sum  of  the  test  times  for  all  of  the 

functions  tested  at  the  particular  echelon.  If  the  test  times 

for  each  function,  T , have  subscripts  which  are  ordered 

m 

according  to  increasing 

c„  a. 

the  test  time  for  the  first  echelon  test  will  be  given  by 
Tv  Wit  ^ Tj, 


wliere  m^  is  the  highest  cost  function  tested  by  the  first 
echelon  and  Tgj^  is  the  time  necessary  for  getting  the  item 
on  and  off  the  tester.  The  second  echelon  test  time  will  be 

c 
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and  the  third  echelon  test  time  will  be 

TL, -p  X, 

where  m2,  m^ , Tq2  and  Tq^  have  meanings  corresponding  to  the 
meanings  of  m^  and  Tqj^.  The  indicies  follow  the  same  rules  for 
second  and  third  echelon  tests  as  for  the  test  times. 


If  the  A are  the  failure  rates  of  the  part  of  the  tester 
m ^ 

that  tests  the  mth  function,  the  failure  rate  of  the  first 
echelon  tester  will  be 


(nT:BF  - 


I 


the  failure  rate  of  the  second  echelon  tester  will  be 


(htbf  ^ 


1*1  L. 


and  the  failure  rate  of  the  third  echelon  tester  will  be 
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where  Aq  is  the  failure  rate  of  the  basic  test  equipment  com- 
prising blowers,  cables,  racks,  and  any  other  equipment  which 
is  shared  by  the  rest  of  the  tester  equipment. 


If  the  A are  the  failure  rates  of  the  part  of  the  tester 
which  tests  the  mth  function  and  the  are  the  repair  times 
for  each  of  these  parts,  the  average  time  to  repair  the 
first  echelon  tester  will  be 

»^t 
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0^  f 
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th.e  average  time  to  repair  the  second  echelon  tester  ^ill  be 

►Vlx. 


•2it 


and  the  average  time  to  repair  the  third  echeloi/ tester 
will  be 


tpA-tr-Tvi^v  t Tt«)f  -s 


where  A is  the  failure  rate  of  the  basic  test  equipment 
0 

and  Rq  is  its  repair  time. 
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COMPUTER  PROGRAM 


The  program  consists  of  a main, program  and  five  subroutines. 
The  purpose  of  each  is: 

MAIN  - controls  program  and  reads  input  data 
INPUTS  - prints  out  the  description  and  value  of  each  data 
element . 

COSTS  - calculates,  for  each  echelon,  the  costs  of  the  test 
equipment  and  of  performing  a test  on  the  item 
CALC  - calculates  the  total  cost  of  testing  the  item  over  its 
service  life. 

XRTNI  - solves  equations  used  in  subroutine  CALC.  XRTNI 
solves  equations  by  using  the  Newton-Raphson 
iterative  technique. 

FCT  - calculates  values  used  by  subroutine  XRTNI 
Table  1 gives  the  computer  program.  The  program  is  in  IBM 
FORTRAN  IV  (Extended)  . 

Operation  of  the  Model 

Figure  3 is  a flow  chart  of  the  program.  Testing  at 
either  two  or  three  echelons  of  test  can  be  simulated.  The 
program  consists  basically  of  four  lo.ops  which  are  executed 
repeatedly.  The  first  two  loops  allow  the  number  of  item 
functions  tested  by  the  first  and  second  echelon  testers  to 
vary  over  their  prescribed  range.  The  remaining  two  loops 
allow  the  intervals  between  tests  at  the  first  and  second 
echelons  to  vary  over  their  prescribed  range. 
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Intervals  between  tests  are  measured  in  use  cycles  which  is  the 
number  of  times  that  the  end  item  has  been  "used'.'.  For  instance, 
a mission  in  which  an  aircraft  takes  off  and  lands  with  a 
particular  missile  constitutes  one  use  cycle  of  that  missile. 

The  purpose  of  subroutine  COSTS  is  to  calculate  the  cost 
of  testing  at  the  ith  echelon  (i=l,2,  or  3)  and  to  calculate 
the  cost  of  the  test  equipment.  Each  cost  depends  only  on 
the  number  of  functions  tested  at  that  level.  For  example, 
suppose  15  functions  are  to  be  analyzed:  the  first  5 at 
echelon  1,  the  next  6 at  echelon  2,  and  the  remaining  4 
at  echelon  3.  Subroutine  COSTS  would  assign  the  cost  of  testing 
functions  1 through  5 to  echelon  1,  of  testing  function^  6 through 
11  ‘to  echelon  2,  and  of  testing  functions  11  through  15  to 
echelon  3.  It  does  this  by  calculating  the  costs  associated 
with  testing  at  each  of  the  three  levels.  Note  that  the  cost 
of  testing  a function  is  charged  to  one  level  only. 

Subroutine  CALC  calculates  the  number  of  use  cycles  that 
elapse  between  tests  using  EQ.2.  This  equation  contains 
nj^  (number  of  cycles  between  first  echelon  tests)  , n2  (number  of 
nj^  cycles  between  second  echelon  tests)  , and  (Number  of  n2 
cycles  between  third  echelon  tests),  as  variables.  If  only 
2 testers  are  used  the  equation  contains  only  nj^  and  n2. 
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For  example  if  n^^lO,  ^2=5,  and  first  echelon  tests  are 

performed  every  n^^^lO  cycles,  second  echelon  tests  every  nj^Xn2= 
10X5=50  cycles,  and  third  echelon  tests  every  nj^Xn2Xn3  = 
10X5X2=100  cycles. 

Due  to  the  complexity  of  the  expression  (Eq.2)  it  is 
solved  for  one  unknown  value  at  a time.  That  is,  given  n^^ 
and  n2  (or  n^  only)  , CALC  will  solve  the  expression  for  n^ 

(or  n2) . The  equation  is  solved  by  using  the  Newton-Raphson 
method,  which  makes  successively  better  approximations.  Given  a 
rough  approximation  to  the  true^  solution  of  the  equation  the 
Newton-Raphson  method,  after  a number  of  iterations,  finds 
a solution  which  is  arbitrarily  close  to  the  true  solution. 

This  procedure  is  carried  out  in  routines  XRTNI  and  FCT. 

After  the  value  of  n^  Cn2  for  the  two  echelon  case)  has 
been  calculated,  the  total  cost  of  testing  at  all  echelons 
over  the  service  life  of  the  item  is  calculated  by  CALC.  The 
results  are  printed.  These  results  are  printed  for  all  values 
of  nj^ , n2,  m^ , and  m2  which  lie  in  the  ranges  specified  by  the 
inputs.  It  is  a simple  matter  to  scan  these  results  to  find 
the  combination  which  will  give  the  lowest  cost  for  the  testing. 


33 


TABLE  2 

COMPUTER  PROGRAM 


A o t 0 U 1 
1 H '1 A K I ^ f O 


Name  calc 


SUHROUT  InE  CALC(Y«PLCtCOST*XNTSTi^*rf<CObT,UStCYC*^sr»IE,U. 
XNFUNC*XNFUN1  tXNFUN?«K  ,L) 

PURPOSE:  TO  CAI.CULATE  THE  CObT(CTEbT)  UF  TEbTINu  ALL  Tnt 
FUNCTIONS  OF  A GIVEN  ITEM  AND  1 Hf  TOIaL  CObT(lLCObT)  OF 
TFbTiNG  THE  ITEM  OVER  I T LIFE  CYCLE.  bU^RUUTINE  CaLT  CALLb 
sugpoutine  xktni  anu  lb  called  HT  maI\j. 

CTEST-COST  OF  TESTING  ITEM  ONE  T 1 Mt 

FNOEF-FF(aCT  ION  OF  ITEMS  IN  USE  wHlCH  '.RE  NOT  DFFECTIVF 
TlCOST-ToTal  COST  OF  TESTING  ITEM  OVER  ITEM  LIFE  CYCLF 
XIEND-MAXIMUM  number  of  ITFRATIUVS  NFinTon-RAPhSoN  METHOD 
allowed  to  GO  THROUGH 

xihtbf-item  Failure  rate 

XNTSI R ( 1 ) -NUMWER  OF  TESTERS  AT  ITH  ECHELON 
XST-INITIAl.  value  USEU  HY  nE  W TUN-R  A RHSO  V METHOD 

COmmON/X/R ( ^ ) 

dimension  TRCOST ( T) tXNTbTR(3) fCOST ( ^) 

equivalence  (XIMTHF*H(1)  ) » (XN1*H(E)  ) * (XNI'*h(3)  ) , (FUEFF1«h(A)  ) . 

« (FOEFFR . 3 (5)  ) » (FuEFF3*B  ( b)  ) . (F  NOEF  « H ( 7 ) ) , ( LN I T EM , H ( « ) ) 
external  FCT 
IFl  AG  = ? 

AFMNI  = < 

AFUN?  = L 

IF(XNFUN?  .NE,  XNF'JNC)  go  to  a 

INTERMEDI ATE  VALUES-TWU  TESTER  CaSE 
Z = XImTBF  » XNl 
C = 1.0 
P = 1.0 

FDEFFJ  = FOEFFV 
GO  ro  S 

intermediate  VALUES-THREE  tester  CASE 
A 7 = XImTsF*  XNl  * XN<? 

C = l.n  - EXP(Z  * (-1.0)  * FOEFFl) 

E = 1.0  - EXP(Z  * (-1.0)  * FDEFF?) 
b F = 1.0  - ExP(XImTHF*  (-1.0)) 

D = 1.0  - EXP(x1mTBF*  (-1.0)  « XNl  FOFFFl) 

const  must  be  >=  I*  OTHERWISE  THE  VALUE  RETURNED  hY  mOUTInE 
XRtnI  will  be  negative.  negative  XNa  OR  XNi  VALUtS  have  no 
meaningful  interpretation. 

CONST  = FnDEF  a (XNl/  FDEFF3)  * F * C / (0  » E) 

IF(CUNST  .GE.  1.0)  GO  TO  EO 
IF(XNFUNE  .EU.  XNFUNC)  go  to  El 
TEMPI  = C0ST(3) 

C0ST(3)  = 0.0 
TEMPE  = THC0ST(3) 

THCOSTO)  = 0.0 
IFLAG  = 1 
GO  TO  35 

El  TEMPI  = COST(E) 

CDST(E)  =0.0  •• 
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NAME  CALC 


LISTING  OF  ESA. LOP. MISSILE. SMCbl 
(CONI ) 


KAOC  00/ 

1 -i  »i  A -•<  I / ft 


r 

c 


r 

r 


r 

r 


o'  » 'f  A * 
NEEO  NEiV 


TEnh?  = TrtCOST(<e) 

TKCuSK?)  = 0.0 
IFl  AG  = 1 
GO  ro  TS 

?n  FPs  = .001 

CALL  XWTNl(X*F.UERF»FCT»XST*EPS»It'iu.it:P.ZliCONST*Z^» 

* XNFUNC.XNFUN?) 

HHANCH  to  APPKOPHIATE  message  IF  NO  V.»LUF  KtTOKNEO  oY  aP'nI 
PS  GO  TO  (100.200) tlEP 

IF(Z1  .LT.  174.  .and.  Z2  .LT.  W4.)  GO  TO  JS 
IF(AnFUN?  .EO.  XnFUNC)  go  TO  2b 
WRITE (b*33)  AFUN1,AFUN2.Xn1*XN2,X 

33  format  ( / .2X  »F  4 .0.1X»F4.0*lX*F4.U.U,F4.0,lX»' 

"'INITIAL  X VALUE  WILL  CAUSE  EXPONENT  UVERFlO>« 

"VALUE*  fSXtFlO. I ) 

GO  TO  41 

2b  WRITE  (6»34)  AFUNl ♦ AFUN2, XNl , X 

34  FORMAT!/. 2X, F4. 0, 1X*F'4. 0»1X. F4. OtiX,  ' " * . 9 X » ' I M I T I AL  X 

"L  Cause  exponent  OVERi-LOw.  need  nEw  starting  VALUE. '.SX 

41  RETURN 

36  T = XIMTRF  " FOEFFl  " XNl 

IF(XmFON?  .NE.  XNFUNC)  GO  (0  4o 


intermediate  QUANTITIES-TwU  tester  case 
S = 0.0 

TRCOSTO)  = 0.(' 

FDEFF3  = 0.0 
X r 0.0 

GO  TO  (44.4  n . IFL  AG 

u=  n.n 
*■  = 0.0 
'-.0  TO  4b 

U = XImtrF  " FOEFFP  " XNl  " XN2 
R = EXR(  (-1.0)  " T " Xn2) 

GO  To  4b 


starting 


VALUE 
. F I 0 . 1 ) 
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46 


INTERMEDIATE 
U = XImTbF  " 


quantities-three  tester  case 

FDEFF2  * XNl  " AN2 


4b 

"3 

4b 


r 

C 


W = FXP(  (-1 .0)  * T " XN2) 

GO  TO  (4b.43) . IFlaG 
S = 0.0 
GO  TO  48 

S r F.XP  ( (-1 .0)  * U " X) 

V = l.n  - E>P( (-1.0)  * T) 

Z = 1.0  - FXP((-1.0)  " U) 

Y = XNITEM  / XM 

calculate  COST  OF  COMPLETELY  TESTI.<lG  ITEM  ONE  TImE 
CTEST  = Y*(COST(l)  ♦ C0ST(2)  " V"  W /(1.0-w))  ♦ (COST(J)/(l 

» *S*V*Y*Z/  (1,0  - S) 

CALCULATE  COST  OF  TESTING  ITEM  OVER  ITS  l.IFECYCLE 
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r"Al7L  aiJi 

LISriNb  OF  tSA. LOR. MISSILE. SKCLl  Irt  HaK  ] S /6 

MEMBER  name  calc  (CUNT) 


TLCOST  = CTEST  YRPlC  « USECYC  ♦ A 1 1 S F R ( 1)  * T RC OS  T ( 1 ) ♦ 

«ANTSTH(P)  * THCUST(P)  ♦ XNTSTR(.3)  **  TR(:uST(i) 

IF(AnFUN?  .EU.  XNFUNC)  bO  TO  70 
IFdFLAG  .NE.  i)  GO  TO  6U 
bO  TO  BO 

BO  WRI  Tt  (^»6b)  AFUNl  »AFUN2»  Aim  i »AN2,FUhFFl  .FOtFF^.FuEFFB,  COST  ( i ) * 

* COST  (?)  .COST (J)  .CTEST, TLCOST 

bb  FORM  T(/.2X»FA.O»1X.F4.0»1A.F4.0»1A»F'*.0.1a»'  1 ’.AX. 

3(F4.?.AX),lX,J(Fa.Q,AX),FO.O.bX,F10.()) 

COST(3)  = TEMPI 
TRCOST(3)  = TEMP? 

return 

n)  IFdFLAG  ,NE.  ?)  GO  TO  7S 
f'-U  TO  BO 

7b  wR  [ TF  (b.  77)  AFUtNl  , aFUIM?,Xim1  ,X  .FOEFH  ,F  DEFF  ? , FOEFF  3,CUST  d ) . 

* COST (2) .COST ( J) »CTEST, TLCOST 

77  format t/,?X.FA,0, 1X.F4.0. lX.FA.O,  lA,  • I « , 1 A . F A . 0 , A A , 

3(F<f.2.<fX)  ,lX.3{Fd.0.‘»A)  ,FR.0,bX,F  l0.()) 
rOST(2)  = TEMPI 
TR0ST(2)  = TEMP? 

RETURN 

BO  WRITE(6.50)  AFUNl,  AF  UN2  , X(M  1 , XN?  , X , FOEF  F 1 , FOLF  F2  , FUEFFU  , COST  (1  ) . 

» COST  (?)  .COST (3)  , CTEST. TLCOST 
SO  format  {/.2X,F4.0.1A,F'.,0.1X,FA.(),1x,F'..(J,1x.F‘.,0,AX, 

* 3(F4,?,AX)  ,1X,3(FB,0.AX)  .FO.O.bX.F  10. 0) 

RETURN 

100  IFdXNF'UN?  .EO.  XNFUNC)  GO  TO  dO 

101  WRTTE(b.l02)  AFUNl .AFUN?.XN1 ,XN2. lEND 

102  FOR -IaT  (/.2X,F4 .0,  lX.FA  .0  . lX.FA.O  , 1 a,F4. 0 , 1 X,  • <>',4a, 

» 'NO  convergence  RITBIN  ',12.*  ITERATIONS*) 

RE  TURN 

110  rtRITE(b»lll)  AFUNl ,AFUN2.XN1 , lENO 
1 1 1 format  (/,?X.FA.0,1X.FA.0»1X,F4.0,1x,i  <»*,9X. 

*>  *NO  convergence  within  *,12, • ITERATIONS') 

RETURN 

200  IF(XnFUN?  .EQ.  XnFUNC)  go  to  210 

201  WRITE(6*202)  AFUNl ,AFUN2,XNl ,XN2 

202  format (/,2X ,F4 . 0 , 1 X»F4 . 0 . IX ,F4. 0 , 1 X ,F4 . 0 , 1 A , • **.4A» 

*»  'AT  SOME  ITERATION  STEP  DERIVATIVE  IS  ZERo*) 

RF  T URN 

210  WRITF(b.211)  AFUn1.AFUN2.Xn1 

211  FORMAT!/. 2X,F4.0,1X.F4,0, lx. F4, 0, lA,  • tt*,9X. 

* *aT  some  iteration  step  OERIVATIVL  is  ZERO') 

RETURN 

ENO 
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wfTMHtW  COSTS 


LISTINb  OF  tSA.LOR.MISSILt.SKCKl 


rAOL  00a 
m "Iak  l-ilh 


SUhPOUT I\E  COSTS (k.L.FOEFFI ♦ FOb FF d * F OE p F T » COS f , AMFUnc » 

* WP  "aTh},  Tt^COST) 

PURPOSE:  TO  CaLCULAIE  Tht  COST  OF  PtRFORM  1 i\iO  i TEST  OM  Tnt 
ITE-'  AT  THE  ITH  ECHElu'M  ( I = 1 »<?*0r<J)  . THE  lEST  COST  aT  A 
ECHEI.OM  depends  On  The  NO'-ihP.w  Ob  FONCPlOUS  OF  THE  I TE 
IHAT  are  CHECKEU  at  TmaT  ECHLLOn. 

CITEM-COST  of  ITE'I 

C0ST(I)-C0ST  of  PEkFOKMINO  test  at  ITh  echelon ( 1=1 3) 

CPACK ( I ) -COST  OF  PACKING  TO  SHIP  ITEM  TO  ITH  ECHELON 

CSHIP ( I ) -COS r OF  SHIPPING  ITEM  10  ITH  ECHELON 

ECOSTS  ( I ) -COST  OF  EOOIPMENT  TO  RlPaIR  I I'H  EONCTION 

EFA  ILR  ( I ) -fa  ILURE  RATE  OK  tOOlPMENT  Tn  kE=»aIR  ITH  FoNCTION 

FDEFFI-PER  CtNf  OF  DEFECTIVE  I TtMS  FAILING  FIRST  ECHEi  on  IEST 

fdeff?-per  cent  of  defective  items  Failing  secoiju  echflon  test 

FDEFFI-PER  CC.NT  OF  defective  ITEMS  FAILING  THIRD  EChElon  TEST 

PAYR{I)-PAY  RATE  AT  ITH  ECHELON ( I = I » ? , J ) 

rpmatr-reRa IR  material  Ratio 

RT IMES ( 1 ) -time  REUD  TO  REPAIR  EOOIP.  TESTING  ITH  FUNCTIUM 
TRCOST-COST  of  ITH  ECHELON  TESTtn 

TTIMES  ( I ) -TIME  RED'D  TO  TEST  THt  ITH  FONCTION  OF  ITF.M 
TURnT ( I ) -TURNaROONU  time  for  item  at  ITh  EChELON 
USECYC-NUMHER  of  ose  cycles 
YHPLC- YEAk:5  per  ITEM  LIFECYCLE 

aIFaIRCD-FaILORE  rate  of  ITH  FONCTION  OF  ITEM 
XMEN(I)-MEN  REO'D  to  operate  TTh  echelon  tester 

COMMON/M/A (3*11) *b(21*b) 

DIMENSION  PaYR(3) *XNmEN(3) *CSH1P(3) ♦CPaCk(3) .TURNT (3) *XNTSTk(3) 
dimension  TEST!  (3)  .XREPRT  ( 3 ) * TRCOS  I ( 3 ) * XJK  m T 3 ( 3 ) * TR  •'■ThF  ( i ) . 

*REPRT  (3)  .C0')T  (3)  ,SOM(3) 

DIMENSION  RT IMES (21) .ECOSTS (21) .FFAlLR(?l).XIFAlR(2I)*TTl3tS(2I) 
equivalence  (PAYR ( 1 ) * A ( 1 * 1 ) ) * (XNmEn ( 1 ) * A ( 1 *2)  ) » (CSHib  ( 1 ) . a ( 1 ,3)  ) . 
* (CPACK ( 1 ) * A ( 1 *4)  ) * ( TURNT ( 1 ) ♦ A ( 1 *S) ) * ( » mTSTR ( I ) ♦ A ( 1 *D)  ) * 

* (ClT'^M.Ad.lD).  (YRPLC*a(2*11)  ) 

equivalence  (RT IMES ( 1 ) *H ( 1 ♦ 1 ) ) * (ECOSTS  ( I)  .3 ( 1 .2)  ) * 

* (EF  AILR  ( I ) .H  ( 1 *3)  ) * ( Alb  aIh  ( 1 ) *3  ( 1 * 4)  ) , 

* ( TTIMES ( I) *H ( I .b) ) 

Nl  = XmFUNC 

KK  = K 
1 = ) 

M=1 

N = XNFUNC  ♦ 1,0 

SUM,TESTT*XREPRT.rRCOST*XTRMTri  ACCOMUlATE  TOTALS  FOR  EACH 
ECHELON  tester  OEPENQlNo  ON  THE  NUHPER  OK  FUNCTIONS  TESfEO 
AT  THAT  level,  COST ( I )* 1=1 *2, 3 IS  IHFN  COMPUTED  USING  THESE 

quantities. 

4 0 SUh(1)=0,() 

TESTT(I)  = TTIMES(N) 

XREPRT(I)=  EFAILRTN)  * HTIMES(N) 

TRCOSTd)  = ECOSTS(N) 

XTPMT8(I)=  EFAILH(N) 
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LISIINIJ  OF  tSA  .LUM.MISbll.E.SKrLl 
mFMHFR  COSTS  (tONJ) 


'^5 

DO  SO  .)  = '^.KK 

SUM(I)  = SUM(l)  ♦ EFaILH(O) 

TcSTTd)  = TFSTT(I)  ♦ Tn''3FS(J) 

XRt'RWT(I)=  XHEPRT(I)  ♦ £FAILR(.l)  « 

TKCUST(I)=  TWCOST(I)  + ECOSTb(J) 

0 u u b b U U ii 
0 00b  700 'I 
OOObHOOO 
OUO  69  0 00 
i)  0 0 b 0 0 0 0 

XTWiV|TO(I)=  XTKMT«(I)  ♦ £FAfL-<(J) 

oouol  000 

">0 

CONT  INOfT 

000b?000 

IF  (I  .GT.  1)  GO  TO  Sb 

OOObbOOO 

1 = 1 + ) 

00  ObA  0(H) 

M=J  ♦ KK 

0 1)  0 b b 0 0 0 

KK  = L 

0 0 0 b 6 0 0 U 

GO  TO  4>0 

00067000 

IF(n1  .NF,  L)  go  to  60 

0 0 0 6 3 0 J i) 

.1  = P 

0006900 0 

COST(3)  = 0,0 

0 0 0 7 0 0 0 0 

GO  TO  70 

0 0 0 7 1 0 0 0 

60 

IF  (I  - 3 .E'J.  0)  GO  TO  6b 

00 J 7?0J0 

r = I ♦ I 

0 0 0 7 3 0 0 0 

M = L ♦ 1 

0007400) 

KK  = W1 

00073000 

GO  TO 

OUO  76OO0 

J = G 

0 0 0 7 7 0 0 0 

70 

DO  7b  1=1. J 

0 0 0 7HoO(i 

TRMTHF(l)  = 1.0/  XTKMTH(I) 

0 0 0 7 9 0 0 (' 

REPKT(I)  = XRtPRT(l)  / SUM(I) 

0 0 0 H 0 0 0 0 

COSTd)  = TESTTd)  * PATH  (I)  ^ ^N-tEN(I)  * CSnlt'd}  ♦ CPACK  d ) 

0 0 0 6 1 0 0 0 

♦ TUPNTd)  « ClTEM/YHPLC  + (PEPRTd)  TESTMd  PaYPCI) 

0 0 0 3 ? 0 0 0 

<» 

♦ Tl-STTd)  » RP'^aTK  * TRCOST(l. ))  / TPMTgF(I) 

OOObUoOO 

7S 

CONTINUE 

00064000 

SU‘' 1=0.0 

ooobboon 

suv<e=o.o 

00066000 

SUv.3  = 0.0 

0006  70  00 

DO  6b  I=1.K 

0 0 0 3 6 0 0 0 

SUMl  = SUMl  ♦ XlFAlHd) 

00  1)69000 

>36 

CONTINUE 

00090000 

DU  90  1=1.1 

00091000 

SUM?  = SUM?  ♦ XlFAlHd) 

0009?000 

un 

continue 

0 0 0 9 3 0 J 0 

IF(N1  ,£3.  1.)  GO  TO  100 

0009<f00  0 

DO  9b  I=1.N1 

000  VbOOO 

SUm3  = SUH3  ♦ XlFAlKd) 

000 96000 

9S 

continue 

0 0 0 9 7 0 0 0 

GO  TO  10b 

0)096000 

1 0 0 

SUM 3 = 0,0 

0 0 0 9 9 0 0 (1 

lOb 

FDLFFI  = SUMl  / XIFAlK(N) 

0 0 i 0 0 0 0 0 

FiJFFFi?  = SUM?  / XIFAIK(N) 

0 0 1 0 1 0 0 0 

FDEFF3  = SUM3  / xlFAlK(N) 

0 0 1 0 ? 0 0 0 

WE TURN 

00103000 

En* 

00  1 04000 
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MFMtJE-W  NA.-'L  ^CT 


listing  UF  ESfi.LOK.MlSSILt.SKCLl 


PAijt  JUb 
i«  -IAr<  i-i  ( t 


SU^■lHOUT  InE  FCT  ( T0L»E  »l»EKF  , «C0NST  » /:?»  ANFJnC»  a mF  uN<e) 

c PUKPOsE:  to  calculate  VALUES  USEU  tiY  SU^WOUTINE  AKTNI. 

r The  F.ilUATlONS  k^HICH  akE  used  DEHENU  0\f  The  yU^-iEki?  Uw  i) 

C UP  echelon  levels  at  which  testing  OCCUhiS.  ECT  IS  CALLtU 

C HY  SUBROUTINE  XHTNI. 

c specifically. FCT  evaluates  the  function  anu  ITS  uEHIVaTIVE 

C Sl^CE  hOTH  if  these  values  are  HEUulWtO  FOh  EaCH  ITERATION  UP 

C the  NEwTON-PAPHbON  METHOD  IN  SUBROUTINE  aKTNI 

r 

COmmon/x/a  (0) 

equivalence  (ALPHA, a ( 1 ) ) , ( XNl , a (P)  ) , ( *N^, a ( J)  ) , (FDEP  P 1 , A (A ) ) , 

» (Fl)EFFP,  A (S)  ) . (P'DEFF  J.  A (b)  ) . (F  NUEE  , A ( M ) , ( XNI  TEM.  A (H)  ) 

if(xnfunp  .eq.  xnfunc)  go  to  b 
GO  TO  b 

r 

C TWO  TESTER  CASE-CHECK  SIZE  OF  EXPONENTS  OCCURInG  In  EQUATIONS, 

r IF  X>174.  E<P(X)  WILL  CAUSE  OVERFLOW  ANO  FRROR  MFSSaGE. 

b Z = ALPHA  * XNI 

71  = I « (-1.0)  TOl 

72  = Z1  » (FOEFFI  ♦ FUEFF2) 

GO  TO  7 

C 

C three  tester  CASE-EXPONENT  CHECK  SlHRlAR  TO  IwO  TESTER  CaSF 

6 7 = ALPHA  <>■  XNI  XN2 

71  = Z * (-1.0)  * TOL 

72  = 71  (1.0  ♦ F0EFF2) 

7 IF(/il  .LT.  174.  .AND.  72  .LT.  17a.)  GO  TO  I 0 
RETURN 

1(1  T = FXP(7  « (-1.0)  * TOD 

IF(XNFijN2  .EQ.  XnFONC)  GO  TO  PG 
C 

C evaluate  derivative  of  FUNCTIU'J  f-thrlp:  tester  case 

S = 1.0  - EXP(7  » (-1.0)  « P'UEFPt'  * TOL) 

20  R = 1.0  - EXP(7  (-1.0)  * TOL) 

OEpF  = (7*H»S  - u*7*FrjEFF2*T*«FDEFFP)  / (S<>S) 

GO  TO  25 
r 

C evaluate  derivative  of  FOnCTIO'I  F-twO  TESTER  CASE 
23  R = l.n  - EXP(Z  * (-1.0)  * FuEFFP  « TOL) 

S = 1.0  - EXP(7  (-1.0)  * FOEFFI  » TOL) 

uEtJF  =(S  * 7 * FOEFF?  » T »»  PDEFF2  - w » 7 « Pl)P.FFl<»T  POEFPl) 
« / (S  S) 

r 

c evaluate  FUNCIIOn  f 
2b  F = (R/S)  - CONST 
RE  TORN 
ENn 


0 0 0 0 1 u 0 0 
00007000 
•J  0 0 0 3 0 0 0 
OUOOnOOO 
0 0 0 0 S 0 J I) 
0 0 0 0 6 0 0 (I 
UOUU  /OOO 
00007100 
0 0 0 0 7 2 0 0 
00007300 
OOOOHUOO 
OOOUROOO 
0 0 0 1 0 0 0 0 
0 0 0 1 1 J 0 0 
00012000 
OOOlzObO 
000X7100 
00012200 
0 0 0 1 3 o 0 0 
00014000 
0 0 0 i o 0 0 0 
000 1 6000 
0 0 0 1 6 0 b 0 
UOOlblOU 
0 0 0 1 7 0 0 0 
0 0 0 1 6 0 0 0 
(JOO  1 ROOO 
0 0 0 2 0 u 0 0 
00071000 
O0'i770on 
00073000 
0 0 0 7 3 0b  0 
0 0 0 7 3 I 0 0 
U 0 07a  0 0 0 
00076000 
00026000 
0007  7000 
OU0770b0 
0007  7 1 on 
00O2HO00 
00070000 
00030000 

0 0 0 3 1 ij  0 n 

000 Jl ObO 
00031100 
Go  0 32000 
IJ  0 0 J 3 0 0 0 
((00  34  00  0 
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J 


ME'IMLH  NiiME  INPUTS 


LISTING  OF  tSr. LOP. MISSILE. SHCt) 


ij  U / 

1b  map  l>7fi 


SUPKOUTlNfc  INPUTS  (XNFUNC»FNUEF»PPMflTR*USECYC»  AST  » AIL  -(D) 
r 

c P'JPPOSE;  TO  UESCPIbt  Each  input  data  ELEmlNI 

r 

COmmON/M/A 
WP I TF ( b» 1 0 ) 

10  FOP'MT (////*S0A. • input  UATA  DESCPIPTION* //) 

*<PirE(bf20)((A(I*J)*I  = itJ)*J=l*ll) 

20  FOP -lAT  ( 10A»  • PAY  PATECF/HP)  AT  I ES  1 1 P • . 0 A . i ( F6 . 2 » 1 1 A ) , / , 

« lOXt'NO.  mE'N  To  OPEPATE  TESIER*  . 19a,  J(E0.2»  l 1A) 

* lOX'COST  OF  SHIPPIimO  TO  T ES  TtR  * » 1 b A , 3 ( F H . 2 » 9X ) , / , 

» lOXt'COST  OF  packing  TO  SHIP  TO  TFS I F P • * 9 A , J ( F 8 . 2 , 9 A ) , / , 

* lox, •TUPNAPOUNO  TIME(YPS)  to  test  1 IES 1 1 a, J (Fb.2» 1 1 a ). /, 

* 10X»«NUMbFP  OF  TESTEpS • »27A» 3 (Fb. 2. 1 1 X)  *///* 

«■  10X,*NO.  TESTEPl  functions  ( INI  r lAL  V AlUE  • . 1 0 « » 3 ( F 3 . 0 , 1 A X ) , 

* /. 1 3X* •Final  value . increment )'♦// » 

* lOX.'NO.  TESTEK2  FUNC T 1 OMS ( INI T I AL  V AlUE • . 1 U A . 3 ( F 3 . 0 f 1 4 A ) , 

« /.13A**FINAL  value. INCPEMtNT) * »//. 

« lOX.'NU.  CYCLES  HETwEEN  FlPil  ECHELON  TtSTS'*7X, 

« 3(F3.0.1AX)  ,/.IJA,  • ( INI  riAL  VAL'it, FINAL  V ALijE  » I NCPEMENT  ) • . 

* //.10X»'NO.  AICYCLES  HETwELN  first  LCHELON  TESTS) '.ax. 
3(F3.0»1AX)  ./.13X, 'CYCLES  mETh/EEN  SECOND  FCHELON  TEbTS',/, 

•»  13X*  • (INITIAL  VALUE.FINAL  v ALUt  , I NCPE  MEN  T ) • * // » 

* lOX.'ITEM  COSr.lTEM  L I F EC  YCLF  . NO  . OF  I T E MS  ' » ^A  , F9  . (J  . 1 2X  , 

« F3.0»1AA»FA.0»//) 

riPITF  (b.?2)  ANFUNC.FnDEF  ,PPt-lATP.iiSECYC»Xsr  t XIE  MO 
22  format ( 10X» 'NO.  ITEM  FUNCTIONS  T E S 1 ED •» 20 X , F3 . 0 .// . 

» 1 OX  , • FPACT  I ON  OF  I IE  AS  IN  USt.  muI  DEF  EC  I 1 VE  ' * 8A  . F 4 , 2 . // . 

« lOX. 'REPAIR  MATERIAL  RA T 1 0 • » 24A t FS . 2 » // » 

* lOX.'NO.  USE  CYCLES' .29X.Fb.O,//, 

* lox. 'Initial  valve  for  jEwTon-raphson  method' *6x»f3.u»//i 

* lOX.'MAAlMljM  NO.  ITERATIONS  FOR  N-R  MET  HUD  ' * bA  . F3 . 0 • //// ) 
NNsXNFUNC  ♦ 1.0 

30  wRlTF(b.31>  (H ( I , 1 ) . 1=1 ,NN) 

31  format  (lux.  • TI-^E  REUO  TO  REPAIR  ENUlP.  WhICH  TESTS  ITH  FUNCTION', 

" fi(Tb7,4(F IS.b.AX)/) ) 

•*0  WPlTE(b,4l)  {B(l,2)  »1  = 1,NN) 

41  FOP -lA  r ( 1 OX , • COS  r of  equipment  to  test  ith  functiom', 

* b(Tb7,4(Flb.6,3X)/) ) 

50  -vRITF(b,51)  (H  ( I ,3)  , 1 = 1 , NN) 

51  FOpMaT ( lox, 'FAILURE  RATE  OF  EUUlR.  TO  TESI  IIH  FUNCTION*. 

* 6(T57,4(r I5,6,3X)/) ) 
bO  WRITE(6,6l)  (H(l,4) ,I=i,NN) 

bi  format  (lox, 'Failure  pate  of  ith  fu-ictiun  of  iieh', 

* b ( Tb7,4 (F Ib.b, 3X) /) ) 

70  ii,PI  TE  (b,71  ) (P  ( I ,5)  , 1 = 1 ,NN) 

71  FOKMAT  (10x,'TIME  RLQD  TU  test  ITH  FUNCTION', 

» b ( T5 7,4 (F 15.6 ,3X) /) ) 

RETURN 

end 


0 0 0 J 1 U 0 0 
00002OUO 
000030  0(1 
0 0 0 0 4 0 0 0 

uootjauoo 

OOOObOOO 
000(1  7(j0  0 
0 0 0 0 H 0 U 0 

oouOpouo 
000  1 uooii 
0 0 0 1 1 0 0 0 
uOU 1 2000 
0 0 0 i 3 0 0 0 
0 0 0 1 4 0 0 0 
0 0 0 i 3 0 0 0 
!)  0 0 i b 0 0 0 
00017000 
0 0 0 i 1 0 0 0 
0 0 (J  I 9 0 0 0 
0 0 0 2 U 0 0 0 
00021000 
000220U0 
00023000 
00024000 
00025000 
0 0 0 2 b 0 0 0 

0 0 02  7 0 0 0 

0 0 0 2 (3  0 0 0 
00029000 
0 0 0 3 0 0 0 0 
0 0 0 3 1 0 0 0 
00032000 
000330UO 
000340 00 
00035000 
0003b000 
0 0 0 3 7 U 0 0 
000  38000 
0003-^000 
0 0 0 4 0 0 0 0 
00041000 
00042000 
UU043UOO 
U0U44000 
00043000 
OOU4bOOO 
0 0 0 4 7 0 0 0 
000480UO 
OOU49UOO 
00030000 
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o c 


ME'-IHEty  NAME  MAIN 


LISTING  OK  tSA. LOR. MISSILE. bhiCEl 


RAuk  OiJM 
IH  Mah  \ ^ th 


C PURPOSE:  HEAD  INPUT  OATA  ANO  CONTRUl  PROGRAM,  THt  KjLLO-«T  JG 

C program  CALC'JLATES  the  life  CtCLE  COST  OK  PERIODICALLY 

r TESTIMG  SPECIFIC  MISSILE  FUNCTIONS.  TESTIxG  IS  aSSU'IKD  TO 

C OCCUR  AT  either  2 OR  A ECHELON  LEVELS.  I HE  PROGRAM  ESSENTIALLY 

CONSISTS  OF  FOUR  LOOPS- HE  FIRST  P '.LLOW  I HE  NUMBER  OF 
FUNCTIONS  TESTED  TO  VARY  OVER  A PRtSCRjfiEJ  RAmGK  AND  THt 
C remaining  2 ALLOW/  THE  TIME  INTERVAL  HETnEEn  TE-STS  TO  V-RY  OVER 

C A PRESCRIBED  RANGE. 

r 

C F 1 STEO-STEPSl/F  FOR  XNFUNl 

C F2STEP-STEPSI2E  FOR  XNK'UN? 

r STEPI-XNl  STEPSIZE 

C STEP2-XN2  STEPSIZE 

C XHVALI-FINAL  XMFUNl  VALUE 

C XHVAL2-FINAL  XNFUN2  VALUE 

C XNl-NO.  CYCLES  RETwIEEN  FIRST  kCHEL.UN  TESTS 

c xnif-final  xNi  value 

r X(n2-N0.  XNI  CYCLES  HETw/EEN  SECOND  ECHELON  TESTS 

C XN2F-FINAL  XN2  VALUE 

^ XnFUNI-NO.  function^  tested  my  first  ECHELON  ( liNll  T 1 AL  VALUE) 

r XnFUNP-NO.  functions  TESTEO  HY  StCUivO  ECHELON  ( iNl  T 1 AL) 

r XNFUNC-NO.  of  item  functions  TEsTEU 

c 


U 0 (J  U 1 0 0 IJ 
(j00iJ2J0() 
1)  U U 0 0 0 U 0 
I)  0 1)  'J  A 0 0 (I 
uOUOdUUO 
U 0 0 U 6 0 V U 
tj  0 U I)  / 0 U 0 

uouunuoo 

U 0 U 0 R 0 0 (1 

u 0 u 1 (j  0 0 0 
000 1 1 uOO 
0 U 0 1 0 U 0 

0 u 0 1 i 0 u 0 
OOOlAOOO 
OU.JISOOI) 
OOUlfaOUO 

uuu  i /’vou 

0 0 'J  1 R 0 0 0 
OOOIRUOO 
U002UU0U 
00021000 
00U2200O 
(>  0 0 2 0 0 0 0 
00024000 


COMHON/M/A  (.3*  1 1 ) »H  (21  .S) 

COHMON/X/C ( W ) 

DI  -ENSION  COST  ( .3 ) , XN  TS  TR  ( 3 ) .TRCOST  (3) 

FQUIVAI  ENCE  ( XNFUNl  » A < 1 » 7)  ) * ( X H V AL  1 . A ( 2 » 7 ) ) * (K  1 STEP  » A ( .3  » 7 ) ) » 
(XNFDN2»A(1»S)  ) » (XriVAL2*A  (2»lO  ) f (F2bT6P»A (3»P)  ) » 
(XMl  .A  ( 1 .0)  ) . (X  glK  , A (2»R)  ) . (STEP]  ♦ A ( .3  » 4 ) ) * 

(XN2*  A ( 1 » 10)  ) * ( XN2F,  A (2*10)),(srEP2tA.(3il0)). 
(XNITF.M.A(3»11)).(XNTSTR(1)»A(1»S)),(  YRPLCt  A (2*  ID) 


WR I T E ( 6 » 1 ) 

formatc  list 


OK  INPUT  DATA*/) 


00O2DU00 
0 0 0 2 S 0 0 0 
0 0 0 2 7 0 0 Cj 
uuudauon 
00024000 
oOOSOOOO 
0 0 0 J 1 0 0 0 
0 0 0 J 2 0 0 0 
0OOJ30O0 
000  34000 
OO03SOU0 


1 1 
15 
C 

FORflAT  (T21  »3  (F10.2.  iX)  ) 
format (T21 *Fl0.2) 

read  ANO  WRITE  Input  data 

DO  5 J=l,ll 

RFAn(5.iO)  (A(I.J)»I=1*3) 

WRITt(b»U)  ( A ( I » J)  * 1 = 1 »3) 

0 0 03 GO  0 0 

0 0 0 J 7 0 0 0 

0 00 .350  0 0 
00034000 

0 004  00  0(1 

0 0 0 4 1 0 0 0 

N 

continue 

00042000 

READ  (b.  1 5)  XNF  UNC  *FnOEF  »RPMATR»UStCYC*x«;T  * xik  nIj 

00043000 

WHI  TE  (R»  15)  XMFUNC*FNOE.F«HPMATR.UStCY(;.XST.XlKNl) 

00044000 

IENO  = XIEND 

0U04S000 

NNsXNFUNC  ♦ 1.0 

0004G000 

DO  40  J=l»5 

00047000 

REAO (5* 100»ENO  = 250)  ( H ( I » J ) * I = 1 f NN ) 

oo04tt(jon 

WRITE (6* 100)  (H  ( I ♦ J)  » I = i *NN) 

0 0 0 4 4 0 0 0 

4 0 

CONI  INUt 

00050000 

100 

F0RMAT(T2lt4Fl5.5) 

00051000 

CALL  INPUTS (XNF UNC* FNDEK  t RPm aTR * UStCYC » X ST  » XIENU) 

00052000 

C ( 1 ) = d (NN*4) 

00053000 

C(7)  * FmOEF 

00054000 

C(M)  = XNITFM 

00033000 

41 


LISTlNij  UP  tb/i.L0R.M[S'^ILE.'5NCtl  IH 

wt'IdtR  NAI'IE  -IAIN  (CONT) 

'JFUNl  = XNPUNl 
NFi)N<?  = XNFUNt* 

IVALl  = XHVALI 
IVal;^  = XHVflL? 

MSTEHl  = FISTEf 
MSTFP<?  = F^STFP 
Nl  = Xixfl 
M2  - X J2 
MF  = XNIF 
M2F  = aNRF 
ISTEP1=  STEPl 
ISTEP2  = STEP2 
'I'MI  IE  (F.*  200) 

200  FOkmaT  (/////t5^X,  tTEbT  COST  SU'^MAHT*, 

* //»2X, 'FUNCTIONS' »5X» 'CYCLES' »7X.  •UE(-ECTI\/E  ITE^S  FaHImO', 

* Sx,  ' TES  I 1 ' » 7X  , ' TEST2'  * 7X.  ' lESTi'  *‘3X,  ' COMPLETE  • *6X, 

« 'I.IFE  CYCLE  ' »/*3Xt  'JMl  ' t3Xt  ' M2' t Ja«  ' N1  ' » JA,  ' m2'  ♦Uxt  'NJ'  *bX» 

« ' TFSTl't2X*'  TEST2',2X»'  TES  T .1 ' ♦ SX  t • COS  T ' . 4X  , ' COS  T ' * hX  , ' C(!S  T » » 

* bXf'TESI  COST ' »BXt 'COST'//) 
r 

r LOOP  ON  FUNCTIONS  TESTED  bY  FIRST  ECHtl.ON  TESTEP 

r 

DO  1000  < = NFUN1  * IVALl  fiNSTEPl 
IF(XnFUN7  .NE.  XNFUNC)  go  TU  IbO 
IVaL?  = MFUN2 
r 

r LOOP  ON  FUNCTIONS  TESTED  bY  SECOND  fcCMFLON  TESTER 

r 

ISO  DO  2000  L=NFUN2» IVAL2.NsrEP2 

if(l  .le.  k)  go  ro  2000 

CALL  COSTS (K»L.FDEFFi ,FDEFF2,P DEFFJt COST, XnFJNC*PPMaTK,THCOSI ) 
C(A)  = FJEFFI 
C(S)  = F0EFF2 
C(b)  = F0EFF3 
f, 

C LOOP  ON  CYCLES  bETwEEN  FIRST  ECHELON  TESTS 

r 

DO  bOO  I=N1 ,N1F , ISTEPl 
XNl  = I 
r (2)  = X >ll 

IF(XNFUN2  .NE.  XMFUNC)  GO  TO  220 
m2F  = m2 

r 

C LOOP  ON  CYCLES  bETWEEN  SECOND  ECHELON  TESTS 

C 

220  00  800  JSN2.N2E, ISTtP2 
XN2  = j 
r(3)  = XM2 

CALL  CALC (YpPLC* COST* XNTSrR»TNCOST,USECYC*XST  * I ENO . XNF UNC • 

* XNFUNl .XNFUN?*K*L) 

800  CONTINUE 

500  continue 
2000  continue 

)000  CONTINUE 

230  STOP  t 


r'l\  jt  ()Ub 
•'  n '<  I T b 


OUObb'lUl) 
0 0 .J  D T 0 0 0 
U 0 0 b b 0 0 0 
(J  0 u b b iJ  0 0 
UOUbOoOO 
0 0 0 o 1 0 0 (J 
0u0b2000 
OOOo  J'iOO 
OUUbAOOO 
0 0 0 b 3 0 u 0 
OOUbbO JO 
0 0 0b  Tool) 
OUObHOUO 
1)  0 0b POOD 
OOUTOOOO 
0 0 0 T 1 0 J 0 
000/2000 
J 0 0 / 3 0 0 0 

oou/<+uuo 

; 0 0 / b 0 0 0 
0 0 0 / b U 0 0 

0 JO  / /ooo 
00 J /bOOO 
0 0 J / b 0 0 0 
0 0 0 d u 0 0 (' 
t)  0 0 b 1 0 0 0 

O00a2oo0 

OOJHJOOu 
OOOdAUUO 
U 0 0 d b 0 0 0 
OOOOl,  JOO 
0 0 J d 7 0 U 0 
oooottool) 
OOObDOOO 

ooosoooo 

000  Vi  000 
0 0 0 V 2 0 0 0 
000D3J00 
OOOVAOOO 
OOOVbOOO 
000  Vb JOO 
0 0 0 y / 0 0 0 
oooydJoo 
oouyvjuo 

OUiuOOOO 
0 0 1 J 1 0 0 0 
00102000 
0 0 1 0 3 0 0 0 

0 0 I 0 u 0 0 
0 J 1 ObOUO 
00 1 ObOOO 
0010/000 
OOlOdOOU 
0 0 1 J y 0 0 0 
001 10000 
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PAbt 

LISTING  OF  ESA.LOS.'^lSSlLt.SHCEl  ' 18  '^AK 

MEMrtEM  NAME  f^AIN  (CONT) 


U 1 U 
1'^  7b 


FNb 


'j  u 1 1 1 0 y 0 


n n 


MFMHEH  NivWE  XKTNl 


LlSriixjb  OF  E'sA.LOR.'^ISSn.fc.SKChl 


-'rtOh.  u 1 1 
IH  man  i-y^h 


SUnROUT  iMt  X^T^'1  (XtF  .0E.RF,FCT  »XST*tRS*ILNO,  I E , Z 1 » CONST  » /£■  * 

^ XnFONC»XNFUN<'‘) 

r 

c PUHROSt:  TO  CALCULATE  THE  VALUE  OF  X USINS  THE  NEwTOn- 

Cf  haphson  technique,  when  testing  is  performed  at  Two 

C ECHELONS  ONLY,  X REPRESENTS  T mE  MUmhe*^  OF  XNl  CYCLES 

C mEIWEEN  second  ECHELON  TESTS.  wHF.m  ( hREE  TYPES  OF  TESTERS 

C ARE  USED,X  REPRESENTS  THt  NUMMtR  OF  XN.  CYCLES  HlTwE'EN 

c Third  echeluim  tests,  xrtni  calls  soHwoijTiNE  fct.  thf.ne  are 

C THRFE  CASES  IN  WHICH  AN  X VAl.iE  IS  NOT  RETURNED. 

C 1)  1ER=1.  NO  CONVERGENCE  IN  lEND  ITtRAlION  STEPS. 

C IER=2.  DIVISION  dY  ZERO. 

C 3)  INITIAL  VALUE  CAUSES  OVERFLOW  iN  CALCULATIONS  INVOLVING 

r exponentials. 

c xrtni  calls  Fct  and  is  called  oy  subroutine  c^lc. 

r 

c prepare  iteration 

IF.R  = 0 
X r XST 
roi  = X 

CALL  FCT  (TUL*F,UERF,Z1  ♦ CONST  , Z(P  , XNF  UNC  * XnK  ON  ^ ) 

TOI.F  = 100.  * EPS 

c start  iteration  loop 

DO  d I=1,IEND 
IF (F) 1 ,7, 1 

C EQUATION  IS  NOT  SATISFIED  UY  X 

] IF (DEPF) 2.b,2 

P IF(aBS(DERF)  - 16.0  (-60.0) ) 8,22, ?2 

c iteration  is  possible 

p?  OX  = F / DERF 

X = X - OX 
TOt.  = X 

TF(aBS(T0L)  .GT.  .1)  GO  TO  12 

X = 1.0 

return 

1?  call  FC T ( TOL *F ,0ERF , Z 1 ,C0NST , Z2 , XNF ONC , XNF On2 ) 

C TEST  ON  SATISFACTORY  ACCURACY 

IF(Z1  .LT.  W4.  .and.  Z2  .IT.  17a.)  Go  TO  lu 
RETURN 

10  TOi  = FPS 
A = AdS(X) 

IF(a  - l.)4,4,3 
3 TOL  = TOL  * A 
A 1F(aBS(DX)  - TODS, b, 6 
B IF(A8S<F)  - TOLF)7,7,d 

6 CONTINUE 

END  OF  iteration  LOOP 

NO  CONVERGENCE  AFTER  ItND  ITERATION  STEPS.  ERROR  RETURv. 

IER  = I 

7 RETURN 

C ERROR  RETURN  IN  CASE  OF  ZERO  DIVISOR 

B lER  a P 
RETURN 
ENO 


0000 1 UUo 
OUOOpOOO 
OOOUJOOO 
0 0 0 0 A 0 0 0 

0000 s 000 

U 0 0 0 ft  0 0 0 

UU0O7U00 
0 0 0 0 ^ 0 0 0 
OOOO'XJOO 
UOVlOuOO 
UOOli'JOO 
0 0 0 i 2 0 0 0 
U 0 0 I 3 0 0 0 
0 0 0 I A V 0 0 
U 0 0 1 ft  0 0 0 
UOOibUOO 
0001 7000 
UOOl  jOOU 
OOO  i >00(1 
U 1)020000 

0002iu00 
UU02?00(.) 
U0023000 
0U02A0  0(I 
0 0 0 2 ft  0 U (y 
00o2b0o0 
0002  7000 
UU02HO00 
OU02QOUU 
00030 000 
0 0 0 J I U 0 0 
0 0 0 J 2 0 0 G 
0 0 V J 3 U 0 0 
O003A0O0 
000 JSOUO 
00  0 S60  00 
0 0 0 d ? 0 0 0 
00  03HUO() 
00039000 
00040000 
OOOAlOOO 
OOOA20UO 
00043000 
00044000 
00048000 
00046000 
00047 000 
00048000 
0004V000 
OOOftOOOO 
0 0 0 ft  1 0 0 0 
00082000 
00083000 
00054000 
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UIST  OF  INPUT  OATA 


TABLE  3 

EXAUPLE  RESULTS  1 


7.00 

12.00 

15.50 

2.00 

2.00 

2.00 

10.00 

500.00 

1000.00 

0.0 

10.00 

10.00 

0.01 

0.10 

0.75 

20.00 

2.00 

1.00 

1.00 

10.00 

1.00 

2.00 

11.00 

2.00 

1.00 

5.00 

2.00 

2.00 

IH.OO 

2.00 

000.00 

10.00 

500.00 

12.00 

0.70 


0.15 

20.00 

7.00 

20.00 

1.000000 

1.500000 

1.000000 

1.500000 

1.500000 

1.500000 

1.500000 

1.500000 

2.000000 

2.000000 

2.000000 

2.000000 

1.000000 

20000.000000 

20000.000000 
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89400.000000 
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100000. oooouo 
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0.000010 
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o.ooiiao 
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0.300000 
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0.300000 
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1.000000 

1.000000 

2.000000 

0.200000 

INPUT  DATA  DESCRIPTION 

PAY  rate(s/hr)  at  tester 

7,00 

12.00 

15.50 

NO.  MEN  TO  OPERATE  TESTER 

2,00 

2.00 

2.00 

COST  OF  SHIPPING  TO  TESTER 

10,00 

500.00 

1000.00 

COST  OF  PACKING  TO  SHIP  TO 

TESTER 

0,0 

10.00 

10.00 

Turnaround  time(yrs)  to  test  item 

0,01 

0.10 

0.75 

NUMBER  OP  TtbTERb 

20.00 

2.00 

1.00 

NO.  TESTERl  FUNC T I ONS (I N I T I AL  VALUE 

1. 

10. 

1. 

final  VALUE, INCREMENT) 

NO.  TFSTER2  FUNCTIONS ( INITIAL  VALUE 

2, 

11. 

2. 

FINAL  VALUE, INCREMENT) 


NO.  cycles  PETKEEN  FIRST  ECHELON  TESTS 
(INITIAL  value, FINAL  VALUE , 1 NCREMENT ) 
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I TST  OF  input  data  EXAV^LE  RESULTS 


NO  CONVFRGENCF  within  2n  ITERATIONS 
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